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ABSTRACT 
Considering  the  absorption  of  terrestrial  radiation  by  carbon  dioxide 
and  water  vapor,  coefficients  of  eddy  conductivity  are  computed  for  heights 
in  the  range  4  to  320  meters  from  data  collected  during  the  Great  Plains 
Field  Program  (August  and  September  1953.)    Power-law  profiles  with  height 
are  fitted  to  the  computed  coefficients  and  values  obtained  from  the  fitted 
profiles  correlated  with  a  stability  parameter. 
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I .        Introduction . 

An  important  step  forward  in  the  theory  of  the  turbulent  transfer  oi 
properties  in  the  lower  layers  of  the  atmosphere  was  taken  when  the  con- 
cept of  Austauschkoeffizienten  or  "exchange  coefficients"  was  proposed 
by  Wilhelm  Schmidt  in  1925  .    Since  that  time  several  attempts  have  been 
made  by  various  workers  to  measure  these  coefficients  through  experi- 
mental observations  in  order  to  test  the  theory  of  turbulent  exchange  = 
In  1950  Cooley  /  5__/  made  determinations  of  eddy  conductivity  near  the 
earth's  surface  from  data  gathered  at  Manor,  Texas.    About  the  same 
time  Jehn  and  Gerhardt  /9_/  made  measurements  of  eddy  heat  flux  up  to 
out  120  meters  at  Manor  and  Gila  Bend,  Arizona .    In  addition  to  these 
were  the  observations  of  Cramer  and  Record  /6__/  near  the  earth's  surface 
in  1953,  and  a  small  number  of  observations  due  to  other  workers  in  the 
period  1955-1957.    These  observations  have  been  summarized  by  Priestly 
/1J/,  and  no  attempt  will  be  made  here  to  review  all  of  the  results  . 

The  primary  reason  that  the  number  of  such  observations  has  been 
ill  is  that  the  accurate  determination  of  micrometeorological  parameters 
is  exceedingly  difficult  and  very  expensive.    For  this  reason  the  quality 
md  volume  of  such  data  obtained  in  the  past  have  been  far  below  the 
level  of  desirability  „ 

The  publication  of  the  "Proceedings  of  the  Great  Plains  Turbulence 
Field  Program,  1  August  to  8  September  1953,  O'Neill,  Nebraska"  made 
new  and  better  data  available  for  the  study  of  the  turbulent  transport  of 
properties  in  the  atmosphere's  lowest  layers.    In  order  to  obtain  new 
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evidence  on  the  nature  of  the  turbulent  heat  and  momentum  transports 

and  their  relationships  to  one  another,  it  was  originally  decided  to 

initiate  thesis  projects  at  the  United  States  Naval  Postgraduate  School  „ 

These  projects  to  be  undertaken  were  to  determine  the  coefficients  of 

(1)  eddy  conductivity,  K„ ,  (2)  eddy  viscosity,  KM  ,  and  (3)  the  ratio 

« 

the  two  coefficients ,  KTT/K      from  this  new  micrometeoro logical  data 

H      M 

In  particular  this  present  study  involves  the  determination  of  the  coeffi- 
cient of  eddy  conductivity  and  the  examination  of  its  variation  with 
height  and  stability  from  the  values  determined. 


2.        The  Great  Plains  data. 

The  data  collection  program  for  the  Great  Plains  project  was  hro'r  en 
iown  into  seven  time-observation  periods  with  a  multitude  of  parameters 
being  observed  at  nearly  equally  spaced  observation  times  within  each 

period . 

For  this  study  the  seventh  observation  period  (7  and  8  September) 

3  used  and  was  considered  to  be  the  most  suitable  of  the  periods  on 
the  basis  of  a  combination  of  two  factors;  (1)  completeness  of  data  uti- 
lized in  the  study  and  (2)  a  minimum  of  advection  during  the  period .    As 

aid  in  the  interpretation  of  the  results  the  following  quotation  /K/ 

is  provided  concerning  the  synoptic  conditions  in  the  area  during  the 

se  enth  period: 

On  7  September,  the  high-pressure  area  east  of  Nebraska  moved 
eastward  while  a  trough  of  low  pressure  approached  from  the  west.    The 
winds  were  southerly  and  initially  light,  becoming  15  to  20  mpho    The 
skies  were  generally  clear  over  Nebraska  with  some  scattered  cirrus 
during  the  morning  of  the  8th.    Cumulonimbus  clouds  and  lightning  were 
visible  to  the  north  of  the  test  site  from  the  afternoon  of  the  7th,  until 
the  early  morning  of  the  8th.    The  seventh  and  last  general  observation 
period  began  at  0405  GST  on  7  September  and  continued  until  1305  GST 
on  8  September  under  relatively  favorable  conditions. 

Of  the  many  tables  of  data  included  in  /10/,  the  following  summary 

shows  those  which  were  used  in  this  study: 

\lejiumber   1  Data  used   2 

Ula  7  September,  1900,  GST,  West  and 

South  components  of  the  geo  atrophic 
wind  at  indicated  pressure  levels. 
Values  estimated  from  local  synoptic 
charts . 


'Table  where  numbers  of  the  form  1.1  etc.  ,  refer  to  tables  in  /10/ 
lere  times  are  omitted  ,  data  were  available  for  each  sounding 


big  Number1  Data  Used   * 

1.3  Station  pressure  and  rela 

2  .  la  Soil  temperature  at  0  -  5  c 

surface  of  the  ground  (tai  en  as  bein  , 

representative  of  the  gro„;.  .  ".-.  r.  . 
ture . ) 

4,2:     p  Mean  air  temperature  fro c    asj 

thermocouples  (MIT,  15-n   nute 
at  2.0,  4.0,  8.0  an  1  I  5 . 0  mete        U  o  « 
the  ground,  UCLA,  average  of 
readings  taken  every  3  minute  3  at 
various  levels  up  to  8.0  metej 
shielded  thermocouples  (JH{  avei 
of  20  readings  at:  each  leve*       ring  a 
5-minute  interval  at  .    el      ip 

to  6.4  meters)  ana  shielded  th€         ;ters 

(UT  ,   10-second  mean..   \~ _••      -■-■    ^h. 

up  to  .8  meters .) 

5.2a  Mean  cross  product  oi   irei  leal  eddy 

velocity  component  and  ill     e    pera- 
ture  multiplied  by  density  ai  .    .  pecific 
heat  of  air  for  3.0,  6,0  and     .         neters 
above  the  ground,    value:  i  ng  for 

1235  and  1435, 

6.4  Radiosonde  pressure,  temperature  and 

relative  humidity  at  sigr.  .    els* 

I  .5  Air  temperature  and  mixii  ilti- 

tudes  from  aerograph  (L-  .      d  ita.   1 
at  start,  top  level  and  end  of  flight  1 
indicated . 

•  3a  Sensible  heat  transfer  to  the  air  from 

the  earth-air  interface  according  to 
theoretical  models  oi  Suon  i  and  Let1  i  . 

During  the  observation  period,  data  were  recorded  bi~hourly  for  even 

hours  of  average  local  time  (CST  plus  35  minutes) .    Sufficient  til  I 

it  i  were  available  to  make  K„  computations  for  six  soundings  0835, 

1033,  1235,  1435  and  1635  on  7  September.    Since  no  night-time  aero  ,.  iph 

ervations  were  made,  only  daytime  determinations  were  possible. 
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2.»  1 .    Instrumental  errors  . 

Contributors  to  the  Great  Plains  project  have  given  the  follow 
values  for  instrumental  errors  inherent  in  the  data; 

1 «,    Micrometeorological  absolute  temperature  errors  (Table  4 


a)  0.2°C 
-b)  0 . 1°C  (MIT  and  UCLA) 
e)  0o05°C 


2 
The  L-20  aerograph  configuration    temperature  and  relative 

errors . 


a)  Temperature  error: 

(i)  1°C  absolute,  for  any  level . 

(2)  Relative  error  between  levels,  0,1  C. 

b)  Relative  humidity  error: 

(1)  to  10%  absolute. 

(2)  Relative  error  between  levels,  3%. 

In  29  cases  a  comparison  of  MIT  temperature  values  at  16  „U  meters 
to  the  L-2Q  data  at  17  meters  gave  a  mean  difference,  MIT  minus  L-20  of 
-  0.10  +  0.11°C. 

For  a  complete  picture  of  the  instrumentation  ,  analysis  of  instru- 
mental errors  and  a  description  of  the  test-site  the  reader  is  referred  to 


"""-Tor  an  interesting  comparison  of  observed  micrometeorological  tempera- 
tures see  fig  .  1  „ 

"The  L-20  configuration  utilizes  the  KS-4  Kollsman  aerograph,  which  is 
an  automatic  device  for  the  automatic  recording  of  temperature  and  rela- 

'    ve  humidity. 
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Fig.  1.  Scatter  of  micrometeorological 
temperature  observations  for  the  0635  sounding. 


Background,  methods  and  equations. 

The  coefficient  of  eddy  conductivity,,  K^,.,  is  defined  by: 

Q 

Kh  - 


where  Q  is  the  turbulent  heat  transfer  in  the  vertical.    It  may  be  observed 
from  this  definition  that  K     is  always  a  positive  quantity.    When  i§       is 
negative  (unstable  conditions) ,  Q  will  be  positive  and  directed  upward „ 
Conversely  for  stable  conditions (st®  >£>)  »  Q  w^  &©  negative  and  directed 
downward . 

In  this  study  starting  values  of  Ktt  at  4  and  8  meters  are  computed 
by  using  eq.  (1)  with  Q  given  by  the  average  value  in  the  layer  £-12 
meters .    For  values  above  the  starting  level  the  following  method  is  utilized 

The  horizontal  turbulent  exchange  of  heat  is  usually  very  small 
compared  to  that  in  the  vertical,  and  vertical  advection  is  usually  eon- 
sidered  to  be  negligible  in  comparison  to  horizontal  advection .    The  fact 
that  clouds  did  not  occur  in  the  layer  under  study  rules  out  a  latent  heat 
contribution ,  and  the  small  height -variation  of  density  can  be  overlooked 
for  the  layers  to  be  considered. 

Then  the  equation  of  turbulent  diffusion  for  heat  becomes! 

Also 


(ty*)™.-  *%t-L**M) 


i 

Close  to  the  ground  the  following  approximations  are  valid: 


t  at     ^  ^i  2V2;  y 


'Hereinafter  the  bars  over  ©  and  T  will  be  omitted  but  understood. 
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)be  io  j, ■    . 

■  :  [i  Jnsliid  tu3         si 

-  u  taoq  ■  ;   ■ 


bojoav  :....:■  '•  too  sJdfiJa  "io:  •■;  . 

> 

: 
b9i.  •  '    "  fbuia 

sifli  e     rl  --  ■  riliw  I   I 

bssil'i;     -  •-  isvsl  LfifiV  10I     .  .• :     • 

'■■■..  r  ■ 

-fii  •  ■      iv  bni  ..         ■      ..''.  oi  ■ 

•   9  ■  ••  >J        ;buia  lebn  .  :>o  ion  :     .  ■  ■•.•''. 

b93lc  dIisv  •-;•'•, 

I  .  ■ 
■   .-  -  srfi  nsii 


bilsv  ,  xoiqqfi  -   txwollol  em  bm..>    | 


•  ■  I  •        ■   '■-        ■■    'o     ted  -j.' . 


Now  returning  to  eq.  (2)  and  substituting  eq.  (3)  we  get; 


fJJ 


5-*T*(&U  +  !iCMS**)] 


(4) 


tr  e     ifference  form  this  becomes: 


■ 


A"*, 


V-V„T-f/^I) 


i  spo  5ing  terms  and  solving  for  K?  gives: 


■ 


x  ^ 


(6) 


3  =  I      The  water  vapor  problem . 

The  water  vapor  radiation  problem  was  approached  using  a  method 
presented  by  D.  L.  Brooks  /!_/  and  based  upon  an  adaptation  of  the 

flux  divergence  equation  of  Bruinenberg: 

r       fit 


Cp    K  loco)       i     £*,     K         n  ^  3*"" 


where  q  is  the  specific  humidity  at  the  reference  level  at  which   — . 

,  to  be  computed;  (rr,.  J    is  the  Elsasser  pressure  correction  factor; 

i 

/is  the  change  in  the  Stefan-Boltzmann  black  body  flux  between 

ecessive  levels  of  temperature  T  in  the  sounding.    The  summation 

index  n  denotes  successively  higher  levels  in  the  atmosphere  above  the 

reference  level  and  the  index  m  the  layers  below  the  reference  level  „ 

1    e  partial  derivative  M? is  the  slope  of  the  curve  of  emissivity  of  water 


ef  "-,-• 


apor,  plotted  against  corrected  optical  depth  and  the  term  &  '      "r~~.     is 

luated  at  the  top  of  the  sounding. 

The  factors d^> are  evaluated  as  finite  difference  slopes  from 
aw/ 

Brooks'  curve    „    In  addition  to  the  curve  Brooks  /  1/  also  presents 

tabular  values  of  the  mean  slope  of  the  emissivity  curve  for  small  incre- 

ments  of  optical  path-length. 

If  the  convention  in  using  eq.  (7)  is  that  ^5-  is  always  positive 

aw 

n  temperature  decreases  with  height  above  the  reference  level,  cooling 

"Not  shown  here  but  presented  in  /  1/. 
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at  the  reference  level  results  from  the  layer  contributions  above , 

.  ly  when  the  temperature  from  the  ground  decreases  up  to  the  ref- 
erence  level, the  contributions  from  below  yield  warming  at  the  reference 
level o    The  second  term  on  the  right,  representing  radiation  divergence 
from  top  of  the  sounding,  gives  a  cooling  contribution  at  the  reference 

level.    The  actual  conventions  used  in  working  with  the  tables  me  some- 
afferent,  however,  and  were  designed  for  maximum  efficiency  in 
at    chine  computation. 

In  the  application  of  the  method,  the  initial  step  is  the  computation 
of  the  corrected  optical  depth,  w,  given  by 

1         (Pi      a    ( '  P    \  ^    i  (8) 


j;  »ferj.. 


where  pQ  and  p  are  the  layer  boundary  pressures  and  pr  is  a  stanc 

pressure,  I00Q  mb„ 

Eq.  (3)  is  approximated  for  a  given  layer  by  using  mean  values  of 
q  and  p,  and  a  finite  pressure  increment  Ap: 

When  increments  of  optical  path  length,  tiw,  are  calculated  for 
each  layer,  values  of  2Jr$may  be  obtained  by  interpolation  in  the  table, 
using  values  w  and  w  +  A  w  in  entering  the  table . 

When  this  work  was  begun  the  decision  was  made  to  perform  the 
computations  by  hand  methods  rather  than  by  electronic  digital  computer, 
IhJ  3  decision  was  in  part  based  on  the  fact  that  use  of  the  then  available 
computer,  the  National  Cash  Register  Company  Model  102A,  would  have 
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necessitated  lengthy  computation  times  exclusive  of  the  time  required 
for  programming.    This  observation  is  based  on  the  fact  that  there  are 
over  1000  entries  in  Brooks1  table  and  only  about  as  many  storage  cells 
in  the  computer.    If  the  102A  had  been  used  it  would  have  been  necessary 
to  store  the  entire  table  on  magnetic  tape,  which  would  virtually  have 
eliminated  the  computation  time-advantage  of  the  computer,    This  diffi- 
culty may  now  have  been  overcome  with  the  advent  of  the  Control  Data 
Corporation  model  1604  computer.    The  larger  storage  capacity  (over 
32,000  cells)  and  extremely  high  speed  of  this  new  computer  would  lend 
itself  well  to  Brooks'  tabular  method. 
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3.2.    The    :arbon    lioxide  >  adiation  problem . 

For  the  carbon  dioxide  part  of  the  radiation  problem  it  was  decided 
to  employ  a  simple  set  of  working  equations  previously  used  by  Brooks 
/  2_/    in  stratospheric  investigations  of  carbon  dioxide  radiative  transfer, 

The  downward  flux  of  infrared  radiation  at  a  level  in  the  atmosphere 
and  for  a  given  frequency  band  (12.2  -  17.5/w  for  CO  J  is  given  by  the 
Schwarzschild  equation  integrated  over  all  angles: 

Fir fo)  =  n(°y?P  (oywt)  4-  £  Bi  a ik  do) 

where  F^{«/.)is  the  downward  flux  of  diffuse  radiation  at  the  reference 

level,  where  »j=wt  ;  w    is  the  optical  path  length  of  C09  above  the  level; 

F  i  (0)  is  the  flux  down  from  the  top  of  the  atmosphere;  "l,=  (0,w,)  is  the 

trans  mi  ssivity  of  diffuse  radiation  by  a  slab  of  optical  path  w  ;  and 

Bi  =  .  185  d  T     is  the  total  black  body  flux    radiated  by  CO2  at  the  approp- 

riate  mean  temperature  of  the  i      layer. 

The  summation  index  i  runs  through  n  layers  from  the  top  of  the  atmosphere 

to  the  level  in  question.  ATpis  the  increment  to  the  flux  trans  mi  ssivity 

at  the  level  w^  contributed  by  each  layer. 

Similarly  the  upward  flux  F    t     (w,)  is: 

Ft  fa)  =  Ff(We)  l^  (*,*/.)  f  2    Bi  A  :,:  (10) 

where  F  t    (w  )  is  the  upward  flux  at  the  ground;  and   't\-    (w,  ,wJ  is  the 

o  h         1      o 

transmissivity  (assumed  independent  of  temperature)  of  the  slab  of 


' 'A  carbon  dioxide  contribution  equal  to  18.5%  of  the  black  body  radiation 
is  quoted  by  various  authors. 
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optical  thickness  w  -w^  ,  or  from  the  ground  to  the  reference  level.    The 
summation  j=l  through  m  runs  from  the  ground  to  the  level  Wi. 

In  using  these  equations  it  is  assumed  that  each  layer  can  be  app- 
roximated by  an  isothermal  layer  whose  absolute  temperature  is  the  actual 
mean  fourth  root  of  the  temperature  of  the  layer  raised  to  the  4th  power. 
Since  the  soundings  do  not  run  to  the  top  of  the  atmosphere,  F   |/  (0)  was 
computed  at  the  top  of  each  sounding. 

The  net  flux  at  the  reference  level  (w  )  is  then  simply: 

FN  (ty)  »  F f  (vt )  -   Fir  (wt)  (12) 

and  the  flux  divergence,  &?„  ,  is  taken  as  being  the  difference  of  net 
flux  at  two  adjacent  levels,  the  assumption  being  that  the  net  flux 
divergence  is  uniform  between  the  two  layers.    The  hourly-temperature 
change  in  terms  of  the  flux  divergence  in  the  layer  is  then: 

(13) 


js  j.         a  on 


2 

where  flux  is  in  cal/cm  /  min  and  pressure  is  in  millibars. 

For  the  determination  of  flux-transmissivity  7^  /  several  choices 
were  available  as  to  method.    Callender  /4_/gives  formulas  for  ^  based 
a  breakdown  of  the  12.5  -  17.5     band  into  five  smaller  wave  length 
intervals.    Burch  and  Howard /3_/,  and  Howard,  Burch  and  Williams 
/_Q_/   give  two  empirical  formulas  based  on  laboratory  data  for  the  total 
absorption  of  the  12.2  -  17.5  j*    band: 

"Strong"  fit:    fA»Jp»   -Mi-G'S  JLt^w  ■+-  +1JU%  (p+f)  (14) 

"Weak"  fit:      JA»  J»  =     3,  /&  w  '^  (p^\  **  (15) 

13 


7f  f*u  ^  -  J  A>  4» 

It  follows  then  that       A  ~     "~s7/Z        ~       "ii"c  c~ZT-' 


and  ^    =  ^   -    A 

In  (14);(15)  w  is  absorber  concentration  in  atmosphere-centimeters;  p 
is  atmospheric  pressure  in  mm  Hg;  and jS  is  the  partial  pressure  of  CO_ 
(assuming  a  constant  ratio     by  weight,   .00046,  of  the  atmosphere.) 

The  Howard  et  al  expressions  (14)  are  for  the  "strong"  fit  or  rela- 
tively large  values  of  w  and  (15)  for  "weak"  fit  or  relatively  small  values 
or  w.    The  cut-off  value  given  is  20%  absorption,  the  log  expression 
applying  for  A  7>  20%  and  the  exponential  expression  applying  for  -<    20%. 

Both  of  the  methods  discussed  above  have  the  disadvantage  of 
extremely  lengthy  computation-time  for  the  data  under  consideration  unless 
an  electronic  digital  computer  is  used.    Furthermore,  the  empirical  equa- 
tion of  Howard  et  al  has  the  disadvantage  of  a  discontinuity  between  the 
two  equations  near  absorptivities  of  20%.    A  third  and  faster  choice  pre- 
sented itself  in  the  form  of  a  curve  giving  fractional  absorption  A  vs  pres- 
sure x  path  length  (mb-cm-atmos)  given  by  Brooks  /_!_/ ,  and  credited  by 
Brooks  to  Burch  and  Howard  /_  3/    (see  fig.  2.)    This  curve  had  the  advan- 
tage of  affording  greater  facility  in  conjunction  with  hand  computations , 
and  it  had  no  obvious  discontinuities.    A  difficulty  presented  itself, 
however,  after  the  C02  computations  were  well  underway  when  it  was 
discovered  that  use  of  the  Brooks  curve  yielded  positive  CO„  radiative 
flux  divergences  near  the  ground  where  large  negative  water- vapor  flux 
divergences  were  obtained.    Since  this  occurred  near  the  level  where 
absorptivities  of  20%  might  be  expected,  the  difficulty  was  overcome  by 
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changing  from  the  curve  to  the  Howard  et  al  "weak"  fit  expression  for 
curve  absorptivities  less  than  20%. 

Although  Brooks  credits  the  curve  to  Burch  and  Howard  /  3/,  it  is 
not  entirely  clear  how  he  obtained  his  curve.    The  Burch  and  Howard  curve 
is  similar  to  that  of  Brooks  except  that  their  ordinate  is  percent  trans- 
mission and  the  abscissa  parameter  is  absorber  concentration  (cm-atmos) . 
The  Burch  and  Howard  curve  was  obtained  from  data  based  on  the  use  of 
the  two  empirical  formulas,  the  discontinuity  being  smoothed  out  in  the 
construction  of  the  curve.    Thus,  regardless  of  how  the  smoothing  was 
done  in  obtaining  Brooks'  curve  it  is  clear  that  it  is  based  on  the  empirical 
expressions  of  Burch  and  Howard.    On  this  basis,  changing  from  the  curve, 
fig.  2,  to  the  "weak"  fit  formula  where  necessary  is  justified,  the  discon- 
tinuity being  accepted . 

For  the  sake  of  comparison,  values  of  percent  absorption  due  to  the 
"weak"  fit  formula  are  plotted  in  fig.  2  along  with  Brooks'  curve.    In  the 
latter,  a  value  of  the  "strong"  fit  formula  corresponding  to  an  cur'e  absorp- 
tion of  20%  is  shown.    It  appears  that  the  smoothing  in  the  region  of 
10-20%  absorption  favored  the  "strong"  fit  curve. 

Since  the  fractional  concentration  of  CO_  by  weight  in  the  atmos- 
phere is  nearly  constant  with  height,  the  amount  near  any  pressure  level 
over  a  given  vertical  distance  may  be  expressed  in  terms  of  a  pressure 
difference.    Thus  every  4.2  mb  of  atmospheric  air  is  equivalent  to  1 
cm-atmos  of  carbon  dioxide  (ie ,  a  layer  of  COo  of  thickness  1  cm  if 
reduced  to  S.T.P.) .    Taking  the  diffuse  radiation  through  a  slab  to  be 
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1.67  times  that  through  a  parallel  radiating  column,  the  equivalent 

optical  path  is: 

V-   ^1       (     * dp  =   .  3CH  Ap  (cm-  *t*ics)  (16) 

AM      Jp,         '  ' 

This  is  the  absorber-concentration  in  the  Burch  and  Howard  expression. 
For  use  with  Brooks'  curve,  w  is  linearly  pressure-weighted  giving  w'  as 
follows: 

*,=  f  iV^-3'^^        (i7> 

with  w"  now  in  mb-cm-atmos  and  p  the  mean  pressure  of  the  layer. 
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4.         Computations. 

4.1.    The  water  vapor  radiative  temperature  change . 

The  radiative  temperature  changes  due  to  water  vapor  were  com- 
puted by  use  of  the  method  discussed  in  section  3.1.    Required  are  the 
knowledge  of  ^temperature ,  pressure,  pressure  differences,  and  relative 
humidity  at  the  given  reference  levels ,  the  reference  levels  being  taken 
as  8,  17,  35,  100,  165,  and  240  meters.    Here? conditions  at  16.0  meters 
are  assumed  to  apply  at  17  meters,  so  that  in  some  cases,  in  order  to 
smooth  out  artificial  inversions  in  the  temperature  data  between  16.0 
and  17  meters,  it  was  necessary  to  shift  the  MIT  temperatures  so  that 
the  16.0  meter  MIT  temperature  would  correspond  to  the  17  meter  L-20 
temperature . 

The  pressure,  and  layer  pressure  differences  were  computed  by 
use  of  the  hypsometric  equation  together  with  the  known  temperature 
profile .    Values  of  specific  humidity  were  obtained  from  the  aerograph 
humidity  data,  radiosonde  data,  and  station  humidity.    Unfortunately 
instrumental  failure  left  a  hole  in  the  micrometeorological  humidity  data 
between  the  station  values  and  the  aerograph  soundings .    For  levels 
between  station  values  and  aerograph  soundings,  linear  interpolation 
with  height  was  assumed  valid  in  arriving  at  values  of  q  for  the  reference 
levels.    As  noted,  radiosonde  data  was  used  to  extend  the  data  above 
the  maximum  level  of  L-20  flight.    This  was  done  by  "hooking"  on  to  the 
nearest  radiosonde  reporting  level  after  the  highest  reported  L-20  data 
were  exhausted. 
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The  computed  water  vapor  radiative  temperature  changes  are 
given  in  table  1 .    Note  the  strong  radiative  warming  tendency  in  the 
lower  layers  at  midday,  tending  to  change  to  cooling  in  the  neighborhood 
of  100  m.    This  effect  has  already  been  noted  by  Cooley  /5_/. 

Table  1.    Infrared  radiative  temperature  change  (°C/Hr.) 
due  to  water  vapor  (times  in  CST) . 


(meters) 

0635 

0835 

1035 

1235 

1435 

1635 

8 

-0.13 

1.25 

2.00 

3.05 

2.23 

0.71 

17 

-0.07 

0.12 

0.17 

1.87 

0.92 

0.10 

35 

-0.02 

0.14 

0.20 

0.13 

0.50 

-0.07 

51 

-0.09 

0.09 

0.07 

-0.11 

0.10 

-0.14 

100 

-0.14 

0.03 

0.02 

0.22 

-0.01 

-0.08 

165 

-0.19 

0.09 

-0.03 

0 

-0.01 

-0.01 

240 

-0.17 

-0.20 

-0.05 

-0.03 

-0.06 

-0.05 
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4.2.    The  radiative  temperature  change  due  to  carbon  dioxide. 

The  CO 2  radiative  temperature  changes  were  computed  by  the 
method  of  section  3.2  and  are  shown  in  table  2.    Since  the  mean  black 
body  flux  <$~X  for  layers  was  needed  in  the  CO?  method,  plots  of  Gi 

versus  height  were  made  for  the  layers  from  0  to  16  meters,  and  0.185x 

_4  4 

<5"T  approximated  by  equal-area  averaging.    Above  16  meters,  0.185<s' T 

was  simply  averaged  between  levels  . 

Since  the  method  yields  (^  vgtjd0„in  layers  between  reference  levels, 

reference-level  values  were  interpolated.    On  the  average,  CO    radiative 

temperature  changes  were  found  to  be  approximately  23%  of  the  water 

vapor  radiative  changes  ,  a  result  which  also  confirms  that  of  Cooley  /5_/ 

in  these  layers.    Note  also  the  tendency  for  the  radiative  temperature  change 

to  change  sign  near  240  meters  in  a  manner  similar  to  that  of  water  vapor. 

Table  2  -.    Infrared  radiative  temperature  change 
due  to  carbon  dioxide  (°C/Hr.) 

Reference  Times 

level  (meters)       0635  0835  1035  1235  1435  1635 

0.26  0.49  0.31  0.13 

0.15  0.26  0.25  0.08 

0.07  0.09  0.10  0.01 

0.05  0.06  0.06  0 

0.02  0.02  0.02  0 

0.01  0.01  0.02  0.01 

0.01  0  -0.01  0 
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8 

-0.10 

0.17 

17 

-0.02 

0.06 

35 

0 

0.03 

51 

0 

0.02 

100 

0 

0.02 

165 

0 

0.01 

240 

-0.06 

0 

4.3.    The  local  rate  of  change  of  temperature. 

The  local  rate  of  change  of  temperature  was  obtained  from  the  MIT 
data  at  4.0,  8.0,  and  16.0  meters  and  from  the  aerograph  temperature 
data  at  35,  51,   100,   165,  and  240  meters.    Since  the  observation  times 
were  centered  bi-hourly  at  35  minutes  after  the  even  hour,  "/)twas  taken 
over  a  two  hour  period  and  centered  at  a  time  35  minutes  after  the  odd 
hour.    For  example,  the  local  change  from  0635-0835  CST  was  centered 
at  0735  CST.    The  aerograph  temperature  changes  were  obtained  in  a 
like  fashion  except  that  some  difficulty  was  encountered  in  the  use  of 
this  data.    A  plot  of  L-20  data  with  height  yields  two  curves  since  data 
is  taken  on  both  ascent  and  descent  of  the  aircraft.    As  an  illustration, 
fig.  3  gives  the  L-20  temperatures  for  the  12  35  sounding.    Where  data 
for  the  given  reference  level  were  missing  smooth-curve  interpolation 
was  used.    This  would  give  an  estimated  "%t  for  a  two-hour  period  either 
from  the  ascent  or  descent  curves  or  both.    Where  doubt  existed  as  to  the 
representative  values  of  ^y^t(i.e.  whether  to  use  the  ascent  or  descent 
curves) ,  either  an  average  value  from  both  ascent  and  descent  curves,  or 
a  subjective  choice  of  the  most  reasonable  sounding  was  used.    Since  all 
periods  were  taken  as  two  hour  periods,  an  error  in  ^y^may  possibly 
exist  due  to  a  displacement  in  time  away  from  the  supposed  observation 
times.    This  type  of  error  is  estimated  to  be  no  more  than  5%  however,  and 
should  be  random  in  nature.    Table  3  presents  the  observed  (or  estimated) 
values  of    v  y^  . 
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Fig.  3.   L-20  aerograph  temperature  sounding  at  123$  CST. 
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Table  3.  Average  local  rate  of  change  of  temperature  (°C/Hr) 

Reference  Times 

level  (meters)  0735  0935  1135  1335  1535 

8  3.24  2.27  1.61  1.07  0.12 

17  2.98  2.23  1.55  1.14  0.15 

35  2.08  2.67  1.41  0.80  0.45 

51  1.58  2.68  1.49  0.67  0.35 

100  0.61  2.61  0.87  0.80  0.40 

165  -0.20  2.15  0.93  0.90  0.37 

240  0.15  1.42  1.05  0.97  0.33 
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4.4.    Advective  temperature  changes  . 

Generally  speaking, in  computing  K„  the  assumption  will  be  made 

rl 

that  advection  is  zero.    However,  it  may  be  shown  from  the  quoted  geo- 
strophic  winds  /lO/  that  this  assumption  is  not  quite  valid  for  the  period 
in  question.    Haltiner  and  Martin  /  7/  have  shown  that  the  layer  mean 
geostrophic  advection  is  given  by: 

4  ZJ  pi 

where  V^.  is  the  layer  mean  geostrophic  wind;  Ty  is  the  layer  mean  virtual 
temperature;  and  AA  is  the  cross-hatched  area  of  the  geostrophic  wind 
hodograph  in  fig.  4. 

Now  computing  the  geostrophic  mean  advection  for  the  layer  950-900 
mb  for  0900  CST  on  7  September  we  get  0.24°C/Hr.    Since  the  actual  wind 
component  perpendicular  to  the  mean  isotherms  is  about  17%  of  the  mean 
geostrophic,  the  advection  term  of  0.24°C/Hr.  will  be  reduced  to 
0.041°C/Hr.    Since  the  wind  speeds  increase  in  the  afternoon  it  is  pro- 
bable that  a  somewhat  larger  advection  will  take  place  at  this  time.    Un- 
fortunately no  other  useful  geostrophic  components  are  given  in  the  data 
other  than  for  0900.    This  consideration  of  advection,  therefore,  requiring 
more  extensive  investigation  will  not  be  pursued  further  in  this  paper    . 


'Except  for  the  0935  sounding  where  the  advection  computed  for  0900 
will  be  used  in  the  determination  of  a  K^for  comparison  with  a  K„ 
computed  without  the  advection  consideration. 
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Fig.  U.  Hodograph  of  the  geostrophic  -wind  for  7  September-0900 
(values  are  estimated  for  the  O'Neill  site  from  local  synoptic 
charts.) 
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4.5.    The  coefficient  of  eddy  conductivity. 

Having  computed  I&J]        it  will  be  necessary  to  obtain  /^J  4^7 
Values  of /^r.+.^j)were  obtained  by  adding  0.0098°C/m  to     AT      . 
To  obtain  AT     above  17  meters  ,  plots  of  aerograph  temperature  with 
height  were  made  on  semi-log  paper  (see  fig.  3). 

Similarly  for  values  up  to  17  meters  plots  of  MIT  micrometerological 
temperatures  with  heights  were  made  on  linear  graph  paper  and  a  smooth 
curve  drawn  through  the  points  .    Values  of    41    were  then  obtained  from 
the  applicable  curve  by  measuring  the  slope  of  the  chord  centered  on  the 
level  in  question.    For  vertical  gradients  at  17  meters     and  below,  one 
meter  intervals  were  taken.    For  the  levels  35  meters  to  165  meters,  IO- 
meter intervals  were  used,  while  20-meter  intervals  were  used  above 
165  meters . 

Returning  to  eq.   (6)  and  omitting  advection ,  one  obtains: 


K- 


Because  of  the  fact  that  AT  was  centered  after  the  odd  hours  and 

fl-fc 

the  radiative  temperature  change  was  centered  after  the  even  hours,  it 
was  decided  to  compute  a  Ktt  centered  after  the  odd  hour  by  doing  the 
following: 

(1)  Averaging  ftT  )         for  two  adjacent  after  "the  even  hour" 
soundings^ giving  a    (MTS\        centered  after  "the  odd  hour". 

■'■•The  16  meter  change  was  assumed  to  apply  at  17  meters. 
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(2)  Repeat  step  (1)  again  for  4T+&  to  obtain  (AT  +  V.  )  , 

values  of  which  are  listed  in  table  4. 

The  previous  steps  yield  the  final  recursion  formula 


cr_     [g-&U3*^»*  (£*»*)■  (21, 


(£**! 


where  K-.  ,  K„  are  also  centered  after  "the  odd  hours". 

Table  4.    Average  lapse  rate  of  potential  temperature,  (°C/meter) . 


Reference 

Time 

levels  (meters) 

0735 

0935 

1135 

1335 

1535 

4 

0.0088 

-0.1011 

-0.1301 

-0.1292 

-0.0992 

8 

0.0173 

-0.0277 

-0.0402 

-0.0482 

-0.0287 

17 

0.0218 

-0.0132 

-0.0252 

-0.0202 

-0.0082 

35 

0.0138 

-0.0137 

-0.0047 

0.0072 

-0.0082 

51                          0.0183         -0.0057  -0.0042         -0.0052         -0.0062 

100                          0.0178           0.0023  -0.0042          -0.0032         -0.0042 

165                          0.0208           0.0073  -0.0027            0.0008            0.0003 

240                          0.0270           0.0143  0.0014            0.0003            0.0018 
As  an  illustration  of  the  form  of  this  equation  for  obtaining  Ktj,  at  17 
meters,  K-,y,  the  equation  is: 

-      [(£J..M-(IU,1&W  *■   &^\                (21) 

/\(7   - 


(£**)„ 


with  A."2  =  17  meters  -  4  meters  =13  meters. 
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Now  a  starting  value  of  K     at  four  meters,  with  temperature  changes 

H 

at  eight  meters,  yields  a  value  of  Ktt  at  17  meters.    Next  to  obtain  the  35- 
meter  Ku  a  starting  value  of  K„  at  eight  meters  with  temperature  changes 
■"■J  rl 

at  17  meters,  is  required.    After  this  the  computed  value  of  K^  at  17  meters 

is  used  to  get  a  value  at  51  meters?  and  the  scheme  is  continued  in  leapfrog 

fashion  until  all  values  have  been  obtained. 

As  previously  noted  in  section  1,  table  5. a  of /10/  gives  values  of 

Q  (meal  cm     min     ).      These  values  were  obtained  from  the  correlation  of 

rapid  fluctuations  of  temperature  and  vertical  motion  using  one-second 

samplings  of  fast  response  probes  together  with  the  Reynold's  formula: 
Q    =     -p   CP   "'  T~' 

As  a  general  principle,  following  Priestly  /H/,  the  turbulent  flux 
of  heat  may  be  considered  constant  with  elevation  in  shallow  layers  above 
about  one  meter.    The  reported  values,  table  5. a  of  /10/,  seem  to  sub- 
stantiate this  result  up  to  12  meters.    Now  values  of  Q  given  in  table  5.2a 
of  /l0_/ are  for  3.0,  6.0  and  12.0  meters.    Averaging  these  to  get  represen- 
tative values  of  Q  in  the  layer  3-12  meters  gives  starting  valu«s  of  Kpj  4 

and  Ktj    3  obtained  by  using  eg.  (1) .  • 

ri  ,  6 

One  further  complication  exists,  however,  in  the  1235  and  1435 
soundings  for  which  no  MIT  values  of  Q  are  given.    Here  it  was  decided 
to  estimate  values  of  Q  from  a  comparison  of  the  MIT  values  of  Q  and 
theoretical  model  values  of  Q0  at  the  earth-air  interface.    Two  such 
theoretical  models     were  available  with  values  for  the  1235  and  1435 

1 . 

For  a  discussion  of  these  theoretical  models  see  Chapter  7.3  of /_10/ 
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_  ,      r  t;S-;  one  due  to  Suomi  /L0_/  and  one  due  to  Lettau  /i0_/. 

Considering  the  entire  seventh  observation  period  Lettau 's  values 
of  QQ  were  found  to  correlate  more  highly  with  the  reported  values  of 
turbulent  heat  transfer  Q  than  did  the  Q     values  of  Suomi.    On  this  basis, 
a  ratio  of  the. average  observed  Q  to  the  Lettau  theoretical  value  Q0  for 
the  1035,   1235,  and  1435  observations  for  all  observation  periods  at 
O'Neill,  Nebraska  was  calculated.     1035  values  were  included  to  give 
enough  values  on  which  to  base  this  ratio.    The  ratio  was  then  applied 
to  the  QQ  of  the  Lettau  theory  for  1235  and  1435  so  that  Q^  and  Qg  for 
these  times  could  be  obtained.    Computed  values  of  Ktt  are  given  in 
table  5 . 
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5.         Results. 

5.1.    Ktt  as  a  function  of  height, 
rl 

Considering  the  values  of  Ktt  given  in  table  5  up  to  the  point  where 
there  is  a  sharp  drop  off  in  values  with  respect  to  height    ,  the  computed 
values  of  K„  *show  a  strong  linear  trend  when  plotted  on  log-log  paper. 

To  substantiate  this  linear  trend. correlation  coefficients  of  log 

Ktt  and  log  *2     were  computed  and  regression  lines  of  best  fit  were  obtained 

and  plotted  for  the  data  for  each  sounding  (see  figs.  5-10) . 

Now  a  straight  line  plot  of  K,,  with  Z  on  log-log  paper  yields  K^ 

ri  rl 

as  a  function  of  exponential  height  according  to  the  general  equation: 

Computed  equations  for  K^  as  an  exponential  function  of  height  are 

ri 

presented  in  table  6. 

Table  6 .    Ktt  as  an  exponential  function  of  height. 


Range  of  Heights  on 

Time  (CST) 

Equation 

which  the  equation  is  based 
(meters) 

0735 

Kw  -  auo  a  '  47Sr 

4  -  320 

0935 

—            __    _              4fro 

4  -  320 

09352 

Kh  *  Fill  *'** 

4  -  320 

1135 

4  -  165 

1335 

KH  =  414.*  ^]AZ* 

4  -  165 

•1535 

KH  -    l£-OG  2. 

4  -  100 

■^Negative  values,  and  values  above  165  meters  for  the  1135  sounding, 

were  excluded . 
^Advection  neglected  in  the  computations  except  for  the  second  result  at 

0935,  for  which  the  0900  computed  advection  was  assumed  to  apply  at  0935. 
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It  is  to  be  noted  that  where  advection  was  considered  (09  35  case) 

the  result  was  not  greatly  different  from  that  in  which  this  parameter  was 

neglected.    The  difference,  however,  indicates  that  the  Ktt-  values, 

with  advection  included,  were  slightly  greater  in  the  lower  portion  of  the 

layer,  and  smaller  in  the  upper  portion. 

It  should  also  be  pointed  out  that  in  most  cases  ,  the  values  of  K„ 

ri 

computed  for  elevations  above  165  meters  showed  a  decrease  below  the 

value  at  the  previous  level.    In  several  cases,  the  values  of  K„  were 

ri 


negative.    This  was  most  likely  due  to  observational  errors  in  V^/ji  -I-&4.)  , 
when  ^y^z-  was  very  close  to  the  adiabatic  lapse,  a  condition  which 
occurred  in  the  neighborhood  of  165  meters  for  the  midday  and  afternoon 


runs 
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5.2.    Kt,  as  a  function  of  stability, 
rl 

It  would  have  been  appropriate  to  compare  values  of  KH  and  m 
to  a  stability  parameter  such  as  the  Richardson  number.    Since  appropriate 

wind  data  for  such  a  comparison  were  not  available,  it  was  more  feasible 

-  2 

to  correlate  K     and  m  with  a  significant  part  of  the  numerator    of  the  in- 
H 

tegrated  Richardson  number  A©  .    This  was  done  and  the  results  are  tabu- 
lated below. 

Table  7.    Linear  correlation  coefficients  of  A©  with  K     andA0  with  m 

at  indicated  heights  . 

K        (meters)        A©  computed  Correlation  Correlation 

over  layers  (meters)  coefficients  of  coefficients  of 

KH  with  A0  m  with  A© 

4  2-4  0.26  -0.73 

8.  2-8  -0.77  -0.71 

17  2-17  -0.96  -0.68 

35  2-35  -0.88  -0.84 

51  2-51  -0.76  -0.75 

100  2-100  -0.67  -0.75 

Considering  values  of  A©  versus  ICj  at  heights  from  four  to  100 

meters  for  all  soundings  (sample  size  30)  a  correlation  coefficient  of 

-0  .43  is  obtained.  <• 


Values  obtained  from  equations  given  in  table  6. 
2Here    A©  =  0fz)-  ©fc) 
where  2  is  variable  and  0(2)  is  the  value  of  ©  at  £    =2  meters 
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6.         Conclusions. 

Although  negative  values  of  K     are  somewhat  disconcerting,  other 

H 

investigators  have  been  confronted  with  this  anomaly  at  much  lower  levels 

/6_/    As  mentioned  earlier,  factors  possibly  contributing  to  these  negative 

values  are  chiefly  non-representative  lapse  rates  due  to  errors  in  the 

measurement  of  temperature,  especially  as  the  lapse  rate  of  temperature 

approaches  the  dry  adiabatic.    Still  another  possibility  arises  when  one 

considers  the  theory  of  Priestly  and  Swinbank  concerning  forced  versus 

free  convection  (see  for  example  p.   261  of /7_/). 

From  an  inspection  of  the  computed  power-law  profiles  for  K^  listed 

rl 

in  table  6,  it  appears  that  m  increases  inversely  with  stability,  the  maxi- 
mum value  of  m  occuring  in  the  early  afternoon.    Although  it  is  doubtful 
whether  the  correlations  listed  in  table  6  can  be  proved  statistically  sig- 
nificant on  the  basis  of  such  a  small  sample  size  (five)  ,  nevertheless  it 
is  felt  that  a  strong  physical  significance  is  inherent  in  the  repetative 

nature  of  the  results  ,  and  that  K     has  a  strong  inverse  dependence  on 

H 

stability.    This  statement  is  given  support  by  the  results  of  a  correlation 
of  the  A0  parameter  with  all  values  of  KH  represented  in  table  7  (sample 
size  30)  where  a  statistically  significant     correlation  coefficient  of  -0.43 
is  obtained.     Further  it  is  apparent,  from  an  inspection  of  table  7  that  this 
correlation  (-0.43)  could  be  further  improved  by  the  elimination  from  con- 
sideration of  values  of  K     at  four  meters  where  a  positive  correlation  of 

H 

K  with  A0  is  obtained.     Considering  the  positive  correlation  of  K 

v 

•At  the  five  percent  level  of  belief. 
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with £ft  and  the  negative  correlation  of  KTT  ,.  with  Aft  ,  it  seems  likely 
that  the  method  of  computing  Ktt  at  four  and  eight  meters  may  have  affected 
the  overall  correlation.  It  must  be  remembered  that  the  assumption  made 
in  this  connection  was  that  Qtor&  is  independent  of  elevation  from  one 
to  12  meters.  At  higher  levels,  the  computational  procedure  gave  con- 
siderably larger  K^-values  which  then  varied  consistently  (if  not  "signi- 
ficantly") with  stability. 

In  summary,  daytime-values  of  Ktt  obtained  under  the  conditions  of 
the  Great  Plains  Turbulence  Field  Program  /1.0/,  tend  to  increase  with  in- 
creasing lapse  rate,  and  with  increasing  height  up  to  at  least  100  meters. 
Above  some  such  level  Kj_j  seems  to  decrease.    It  is  possible  that  a 
specific  function   /^  —  -f(Rc)      may  exist,  but  the  data  was  not  sufficiently 
detailed  to  establish  such  a  result  at  this  time. 
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